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Abstract 

One approach to the analysis of stochastic fluctuations in market prices is to model charac- 
teristics of investor behaviour and the complex interactions between market participants, with 
the aim of extracting consequences in the aggregate. This agent-based viewpoint in finance goes 
back at least to the work of Garman (1976) and shares the philosophy of statistical mechanics in 
the physical sciences. We discuss recent developments in market microstructure models. They 
are capable, often through numerical simulations, to explain many stylized facts like the emer- 
gence of herding behavior, volatility clustering and fat tailed returns distributions. They are 
typically queueing-type models, that is, models of order flows, in contrast to classical economic 
equilibrium theories of utility-maximizing, rational, "representative" investors. Mathematically, 
they are analyzed using tools of functional central limit theorems, strong approximations and 
weak convergence. Our main examples focus on investor inertia, a trait that is well-documented, 
among other behavioral qualities, and modelled using semi-Markov switching processes. In par- 
ticular, we show how inertia may lead to the phenomenon of long-range dependence in stock 
prices. 
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1 Introduction 

Modeling market microstructure in order to understand the effects of many individual investors 
on aggregate demand and price formation is both a classical area of study in economics, and a 
rapidly growing activity among researchers from a variety of disciplines, partly due to modern- 
day computational power for large-scale simulations, and the increased availability of price and 
order-book data. Among the benefits of this type of analysis, whether mathematical or simulation- 
based, is the design of better models of macroscopic financial variables such as prices, informed by 
microscopic (investor-level) features, that can then be utilized for improved forecasts, investment 
and policy decisions. 

The approach we discuss here is to identify characteristics common to large groups of investors, 
for example prolonged inactivity or inertia, and study the resulting price dynamics created by order 
flows. Typically, we are interested in understanding the microstructure effects on the aggregate 
quantity through approximations from stochastic process limit theorems when there is a large 
number of investors. 

In this point of view, we model right away the behavior of individual traders rather than char- 
acterizing agents' investment decisions as solutions to individual utility maximization problems. 
Such an approach has also been taken in [45], [40], [65] and [39J, for example. As pointed out by 
O'Hara in her influential book Market Microstructure Theory [72], it was Garman's 1976 paper 
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[45j that "inaugurated the explicit study of market micro-structure" . There, he explains the phi- 
losophy of this approach as follows: "The usual development here would be to start with a theory 
of individual choice. Such a theory would probably include the assumption of a stochastic income 
stream [and] probabilistic budget constraints • • • . But here we are concerned rather with aggregate 
market behavior and shall adopt the attitude of the physicist who cares not whether his individual 
particles possess rationality, free will, blind ignorance or whatever, as long as his statistical me- 
chanics will accurately describe the behavior of large ensembles of those particles" . This approach 
is also common in much of the econophysics literature (see the discussion in [36] , for example) , and 
is of course prevalent in queueing models of telephone calls or internet traffic [20J , where interest is 
not so much on causes of phone calls or bandwidth demand, but on phenomenological models and 
their overall implications. As one econophysicist explained it in reaction to the usual battle-cry of 
the classical economist about rational behaviour, when AT&T uses queuing models, it doesn't ask 
why you call your grandmother. 

In this article, we provide an outline to recent surveys on agent-based computational models 
and analytical models based on dynamical systems, while our focus is on developing limit theorems 
for queueing models of investor behaviour, which apply modern methods from stochastic analysis to 
models based on economic intuition and empirical evidence. The goal is in obtaining insights into 
market dynamics by understanding price formation from typical behavioral qualities of individual 
investors. 

The remainder of this paper is summarized as follows: in Section [2j we briefly survey some 
recent research on agent based models. These models relate the behavioral qualities of investors 
and quantitative features of the stock price process. We give a relevant literature review of Queuing 
Theory approaches to the modeling of stock price dynamics in Section [2T3T In Section [2~4"l we discuss 
evidence of investor inertia in financial markets, and we study its effect on stock price dynamics 
in Section [3l Key tools are a functional central limit theorem for semi-Markov processes and 
approximation results for integrals with respect to fractional Brownian motion, that establish a 
link between investor inertia and long range dependence in stock price returns. These are extended 
in Section 13.21 to allow for the feedback of price of the stock into agents' investment decisions, 
using methods and techniques from state dependent queuing networks. We establish approximation 
results for the stock price in a non-Walrasian framework in which the order rates of the agents 
depend on the stock price and exogenously specified investor sentiment. Section [4] concludes and 
discusses future directions. 

2 Agent-Based Models of Financial Markets 

In mathematical finance, the dynamics of asset prices are usually modelled by trajectories of some 
exogenously specified stochastic process defined on an underlying probability space (f2,.F, P). Geo- 
metric Brownian motion has long become the canonical reference model of financial price dynamics. 
Since prices are generated by the demand of market participants, it is of interest to support such 
an approach by a microeconomic model of interacting agents. In recent years there has been in- 
creasing interest in agent-based models of financial markets where the demand for a risky asset 
comes from many agents with interacting preferences and expectations. These models are capable 
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of reproducing, often through simulations, many "stylized facts" like the emergence of herding 
behavior [39j 65] ; volatility clustering (2U [66] , or fat-tailed distributions of stock returns [25] , that 
are observed in financial data. 

In contrast to the traditional framework of an economy with a utility-maximizing representa- 
tive agent, agent-based models typically comprise many heterogeneous traders who are so-called 
boundedly rational. Behavioral finance models assume that market participants do not necessarily 
share identical expectations about the future evolution of asset prices or assessments about a stock's 
fundamental value. Instead, agents are allowed to use rule of thumb strategies when making their 
investment decisions and to switch randomly between them as time passes. Following the seminal 
article of Frankel and Froot |41j . one typically distinguishes fundamentalists, noise traders and 
chartist^]. A fundamentalist bases his forecasts of future asset prices upon market fundamentals 
and economic factors such as dividends, quarterly earnings or GDP growth rates. He invests in 
assets he considers undervalued, that is, he invests in assets whose price is beneath his subjective 
assessment of the fundamental value. Chartists, on the other hand, base their trading strategy 
upon observed historical price patterns such as trends. Technical traders try to extrapolate future 
asset price movements from past observations. Fundamentalists and chartists typically coexist with 
fractions varying over time as agents are allowed to change their strategies in reaction to either the 
strategies' performances or the choices of other market participants. Some of these changes can be 
self reinforcing when agents tend to follow the more successful strategies or agents. This may lead 
to temporary deviations of prices from the benchmark fundamental or rational expectations prices 
generating bubbles or crashes in periods when technical trading predominates. Fundamentalists 
typically have a stabilizing impact on stock prices. 

In this section, we review some agent-based models of financial markets. Our focus will be on a 
class probabilistic models in which asset price dynamics are modelled as stochastic processes in a 
random environment of investor sentiment. These models are perhaps most amenable to rigorous 
mathematical results. Behavioral finance models based on deterministic dynamical systems are 
covered only briefly as they are discussed extensively in a recent survey by Hommes [49] , For 
results on evolutionary dynamics of financial markets we refer to [47] . [32], or [75] and references 
therein. 

2.1 Stock Prices as Temporary Equilibria in Random Media 

Follmer and Schweizer [40] argue that asset prices should be viewed as a sequence of temporary 
equilibrium prices in a random environment of investor sentiment; see also |37j . In reaction to a 
proposed price p in period t, agent a G A forms a random excess demand e®(p, u), and the actual 
asset price Pt(u>) is determined by the market clearing condition of zero total excess demand. 
In [30], individual excess demand involves some exogenous liquidity demand and an endogenous 
amount obtained by comparing the proposed price p with some reference level Sf. This dependence 
is linear on a logarithmic scale and individual excess demand takes the form 



Survey data showing the importance of chartist trading rules among financial practitioners can be found in, e.g., 




(2.1) 



[81] and [42] 
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with nonnegative random coefficients c"(u;). Here rjf(w) is the individual's liquidity demand. The 
logarithmic equilibrium price St(u)) := log Pt(uj) is then determined via the market clearing condi- 
tion XlaeA e a (i-^ (oj) , oj) = 0. It is thus formed from an aggregate of individual price assessments and 
liquidity demands. If the agents have no sense of the direction of the market and simply take the 
last logarithmic price St—i as their reference level, i.e., if logSf = St-i, then the log-price dynamics 
reduces to an equation of the form 

S t = S t -i + r] t 

were rjt denotes the aggregate liquidity demand. In this case the dynamics of logarithmic prices 
reduces to a simple random walk model if the aggregate liquidity demand is independent and 
identically distributed over time. This is just the discretized version of the Black-Scholes-Samuelson 
geometric Brownian motion model. 

A fundamentalists bases his investment decision on the idea that asset prices will move closer 
to his subjective benchmark fundamental value F a . In a simple log-linear case 

log^-SU + (2.2) 

for some random coefficient < a" < 1. If only such information traders are active on the market, 
the resulting logarithmic stock price process takes the form of a mean-reverting random walk in 
the random environment {at}j g N {o-t = {«"}aeA)- A combination of information trading and 
a simple form of noise trading where some agents take the proposed price seriously as a signal 
about the underlying fundamental value replacing F a in (|2.2p by p leads to a class of discrete time 
Ornstein-Uhlenbeck processes. Assuming for simplicity that subjective fundamentals equal zero 
the logarithmic price process takes the form 

S t - S t -i = nftSt-i + It (2.3) 

with random coefficients jt an d It- These coefficients describe the fluctuations in the proportion 
between fundamentalist and noise traders. When noise trading predominates, % becomes negative 
and the price process transient. Asset prices behave in a stable manner when the majority of the 
agents adopts a fundamentalist benchmark. 

2.1.1 Random environment driven by interactive Markov processes 

Let us now discuss a possible source of randomness driving the environment for the evolution of 
stock prices. Extending an earlier approach in [37], Horst [53] analyzes a situation with countably 
many agents located on some integer lattice A where the environment is generated by an underlying 
Markov chain with an interactive dynamics. There is a set C of possible characteristics or trading 
strategies. An agent's state xf € C specifies her reference level for the following period. The 
environment is then driven by a Markov chain 

U(x t ;-) = n^(x t ;-) (2.4) 

where xt = (x" ) a eA denotes the current configuration of reference levels. The distribution TT a (xt; •) of 
an agent's state in the following period may depend both on the current states of some "neighbors" 
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and signals about the aggregate behavior. Information about aggregate behavior is carried in 
the empirical distribution g{x t ) or, more generally, the empirical field R(xt) associated to the 
configuration x%. The empirical field is defined as the weak limit 



along an increasing sequence of finite sub-populations A n j A and (0 a ) a eA denotes the canonical 
shift group on the space of all configurations. Due to the dependence of the transition probabilities 
Tr a (x; •) on aggregate behavior, the kernel IT does not have the Feller property, and so standard 
convergence results for Feller processes on compact state spaces do not apply. As shown by Follmer 
and Horst [38] and Horst [50J the evolution of aggregate behavior on the level of empirical fields 
can be described by a Markov chain. In [53], it is the {i?(x^)}t g N process that generates the 
environment: 



Under a weak interaction condition the process {72(xt)}tgN settles down to a unique limiting distri- 
bution. Hence asset prices asymptotically evolve in a stationary and ergodic random environment. 
This allows us to approximate the discrete time process {StjteN by the unique strong solution to 
the stochastic differential equation 

dZ t = ZtdXt + dX t 

where X and X are Brownian motions with drift and volatility; see [30] or [53] for details. 
2.1.2 Feedback effects 

The random environment in [53] is generated by a Markov process describing the evolution of 
individual behavior. While this approach is capable of capturing some interaction and imitation 
effects such as word-of-mouth advertising unrelated to market events, the dynamics of {xt}t^ lacks 
a dependence on asset price dynamics. The model by Follmer, Horst, and Kirman [39J captures 
feedback effects from stock prices into the environment. At the same time it allows for trend 
chasing. A trend chaser or chartist bases his expectation of future asset prices and hence his 
trading strategy upon observed historical price patterns such as trends. In [39], for instance, the 
chartist's benchmark level takes the form 



A combination of the trading strategies (|2.2|) and (|2.5|) yields a class of asset price processes that can 
be described by a higher order stochastic difference equation. In [3U], the agents use one of a number 
of predictors which they obtain from financial "gurus" to forecast future price movements. The 
agents evaluate the gurus' performance over time. Performances are measured by weighted sums 
of past profits the strategies generate. The probability of choosing a given guru is related to the 
guru's success. As a result, the configuration x% of individual choices at time t is a random variable 
whose distribution depends on the current vector of performance levels Ut-\- This dependence 




(7t)7t) ~ Z(R(xt); •) for a suitable stochastic kernel Z. 



St-i + Pt(St-i — St-2)- 



(2.5) 
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of the agents' choices on performances introduces a feedback from past prices into the random 
environment. Loosely speaking one obtains a difference equation of the form (|2.3I) where 

(7t>7i) ~ Z{Ut', •) for a suitable stochastic kernel Z. 

While prices can temporarily deviate from fundamental values, the main result in [39] shows that 
the price process has a unique stationary distribution, and time averages converge to their expected 
value under the stationary measure if the impact of trend chasing is weak enough. 

2.1.3 Multiplicity of equilibria 

As argued by Kirman [58], in a random economy with many heterogenous agents, a natural idea of 
an equilibrium is not a particular state, but rather a distribution of states reflecting the proportion 
of time the economy spends in each of the states. In the context of microstructure models where 
liquidity trading or interaction effects prevent asset prices from converging pathwise to some steady 
state, stationary distributions for asset prices are thus a natural notion of equilibrium. In this sense, 
the main result in [39] may be viewed as an existence and uniqueness result for equilibria in financial 
markets with heterogenous agents. Horst and Wenzelburger [53] study a related model with many 
small investors where performances are evaluated according historic returns or Sharpe ratios. In the 
limit of an infinite set of agents the dynamics of asset prices can be described by a path dependent 
linear stochastic difference equation of the form 

y t = i4(ft_i)y t _i + B( ft _i,e t ). 

Here {e^teN is an exogenous i.i.d sequence of noise trader demand and Qt-i denotes the empirical 
distribution of the random vector Yq, Y\, . . . , Yt-\- While the models shares many of the qualitative 
features of [53] and [39], it allows for multiple limiting distributions of asset prices. If the interac- 
tion between different agents is strong enough, asset prices converge in distribution to a random 
limiting measure. Randomness in the limiting distribution may be viewed as a form of market 
incompleteness generated by contagious interaction effects. 

2.1.4 Interacting agent models in an overlapping generations framework 

The work in [54] is based on earlier work by Bohm et al. |12j . Bohm and Wenzelburger [13] . 
and Wenzelburger [82]. These authors developed a dynamic analysis of endogenous asset price 
formations in the context of overlapping generations economies where agents live for two periods 
and the demand for the risky asset comes from young households. They investigate the impact 
of different forecasting rules on both asset price and wealth dynamics under the assumption that 
agents are myopic and therefore boundedly rational, mean-variance maximizers. Bohm et al. [12] 
study asset prices and equity premia for a parameterized class of examples and investigate the role 
of risk aversion and of subjective as well as rational beliefs. It is argued that realistic parameter 
values explain Mehra and Prescott's equity premium puzzle (|70j). The model is generalized in [82] 
to a model with an arbitrary number of risky assets and heterogeneous beliefs, thus generalizing the 
classical CAPM. A major result is conditions under which a learning scheme converges to rational 
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expectations for one investor while other investors have non-rational beliefs. A second major result 
is the notion of a modified market portfolio along with a generalization of the security market line 
result stating that in a world of heterogeneous myopic investors, modified market portfolios are 
mean-variance efficient in the classical sense of CAPM, regardless of the diversity of beliefs of 
other agents. See [TT] for a related approach. 

2.1.5 Feedback Effects from Program Trading, Large Agents and Illiquidity 

A different type of feedback effect, from the actions of a large group of program traders or large 
influential agents has been modelled in the financial mathematics literature. In the 1990s, following 
the Brady report that attributed part of the cause of the 1987 stock market crash to program trading 
by institutions following portfolio insurance strategies, researchers analyzed the feedback effect 
from option Delta-hedging by a significant fraction of market participants on the price dynamics 
of the underlying security. See, for example, Prey and Stremme [35] . Sircar and Papanicolaou [75], 
Schonbucher and Wilmott [76] and Platen and Schweizer [73] . 

Related analyses can be found in models where there is a large investor whose actions move the 
price, for example Jonsson and Keppo [55], and where there is a market depth function describing 
the impact of order size on price, for example Cetin et al. [18J. A cautionary note on all such 
models is that, under sensible conditions, they do not explain the implied volatility smile/skew 
that is observed in modern options markets (in fact they predict a reverse smile). This would 
suggest that program trading, large agent or illiquidity effects are second order phenomena as far 
as derivatives markets are concerned, compared with the impacts of jumps or stochastic volatility. 

There has also been some recent empirical work on estimating the market depth function, in 
particular the tail of the distribution governing how order size impacts trading price: see Farmer 
and Lillo [35] and Gabaix et al. [44| . 

2.2 Stock Prices and Random Dynamical Systems 

An important branch of the literature on agent-based financial market models analyzes financial 
markets in which the dynamics of asset prices can be described by a deterministic dynamical system. 
The idea is to view agent-based models as highly nonlinear deterministic dynamical systems and 
markets as complex adaptive systems, with the evolution of expectations and trading strategies 
coupled to market dynamics. Many such models, when simulated, generate time paths of prices 
which switch from one expectations regime to another generating rational routes to randomness, 
i.e., chaotic price fluctuations. As these models are considerably more complex than the ones 
reviewed in the previous section, analytical characterizations of asset price processes are typically 
not available. However, when simulated, these model generate much more realistic time paths of 
prices explaining many of the stylized facts observed in real financial markets. 

Particularly relevant contributions include the early work of Day and Huang [27] , Frankel and 
Froot [H] and the work of Brock and Hommes p3]. The latter studies a model in which boundedly 
rational agents can use one of two forecasting rules or investment strategies. One of them is costly 
but when all agents use it, the emerging price process is stable. The other is cheaper but when 
used by many individuals induces unstable behavior of the price process. Their model has periods 
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of stability interspersed with bubble like behavior. In [15] the same authors introduced the notion 
of Adaptive Belief Systems (ABS), a "financial market application of the evolutionary selection 
of expectation rules" analyzed in [14j . An ABS may be viewed as asset pricing models derived 
form mean- variance optimization with heterogenous beliefs. As pointed out in [Hj], "a convenient 
feature of an ABS is that it can be formulated in terms of (price) deviations from a benchmark 
fundamental and (...) can therefore be used in experimental and empirical testing of deviations 
from the (rational expectations) benchmark." Recently, several modifications of ABSs have been 
studied. While in [15] the demand for a risky asset comes from agents with constant absolute 
risk aversion utility functions and the number of trader types is small, Chiarella and He [21j and 
Brock, Homines, and Wagener |16] developed models of interaction of portfolio decisions and wealth 
dynamics with heterogeneous agents whose preferences are described by logarithmic CRRA utility 
functions and many types of traders, respectively. Gaunersdorfer [46] extends the work in p~4] to 
the case of time-varying expectations about variances of conditional stock returnj§. 

2.3 Queuing Models and Order Book Dynamics 

The aforementioned models differ considerably in their degree of complexity and analytical tractabil- 
ity, but they are all based on the idea that asset price fluctuations can be described by a sequence 
of temporary price equilibria. All agents submit their demand schedule to a market maker who 
matches individual demands in such a way that markets clear. While such an approach is consis- 
tent with dynamic microeconomic theory, it should only be viewed as a first steps towards a more 
realistic modelling of asset price formation in large financial markets. In real markets, buying and 
selling orders arrive at different points in time, and so the economic paradigm that a Walrasian 
auctioneer can set prices such that the markets clear at the end of each trading period typically 
does not apply. In fact, almost all automated financial trading systems function as continuous 
double auctions. They are based on electronic order books in which all unexecuted limit orders 
are stored and displayed while awaiting execution. While analytically tractable models of order 
book dynamics would be of considerable value, their development has been hindered by the in- 
herent complexity of limit order markets. So far, rigorous mathematical results have only been 
established under rather restrictive assumptions on aggregate order flows by, e.g., Mendelson [71] . 
Luckock [62j and Kruk [60]. Statistical properties of continuous double auctions are often analyzed 
in the econophysics literature e.g., Smith et al. [79] and references therein. 

Microstructure models with asynchronous order arrivals where orders are executed immediately 
rather than awaiting the arrival of a matching order and where asset prices move into the order 
to market imbalance are studied by, e.g. Garman [35]; Lux [631 E3 EI] or Bayraktar et al. ([5]). 
These models may be viewed as an intermediate step towards a more realistic modeling of electronic 
trading systems. 

A convenient mathematical framework for such models, which we will develop in detail in 
Section [372], is based on the theory of state-dependent queuing networks (see [68] or [69] for detailed 
discussions of Markovian queuing networks) . Underlying this approach is the idea that the dynamics 

2 There are many other papers utilizing dynamical system theory to analyze asset price dynamics in behavioral 
finance models. For a detailed survey, we refer the interested reader to [49] . 
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of order arrivals follows a Poisson-type process with price dependent rates and that a buying (selling) 
order increases (decreases) the stock price by a fixed amount (on a possibly logarithmic scale to 
avoid negative prices). 

More precisely, the arrival times of aggregate buying and selling orders are specified by indepen- 
dent Poisson processes II + and n_ with price and time dependent rates A + and A_, respectively, 
that may also depend on investor characteristics or random economic fundamentals. In the simplest 
case the logarithmic price process {St}t>o takes the form 

St = S + ~tt+(^J X + (S u ,u)du\ - II (^j X-(S u ,u)du\ . 

The excess order rate X+(S U , u) — A_ (S u , u) may be viewed as a measure of aggregate excess demand 
while n + (Jq X+(S u , u)du\ — II_ ^ J Q * A_(5 U , u)duj denotes the accumulated net order flow up to 
time t. In a model with many agents and after suitable rescaling the asset price process may be 
approximated by a deterministic process while the fluctuations around this first order approximation 
can typically be described by an Ornstein-Uhlenbeck diffusion. 

Recently, such queuing models have also been applied to modeling the credit risk of large 
portfolios by Davis and Esparragoza [26]. They approximate evolution of the loss distribution 
of a large portfolio of credit instruments over time. We further elaborate on queuing theoretic 
approaches to stock price dynamics in Section [3l Before that, we introduce a common investor 
trait, investor inertia, and show the effects of this common trait on stock prices. 

2.4 Inertia in Financial Markets 

The models mentioned previously assume that agents trade the asset in each period. At the end 
of each trading interval, agents update their expectations for the future evolution of the stock 
price and formulate their excess demand for the following period. However, small investors are 
not so efficient in their investment decisions: they are typically inactive and actually trade only 
occasionally. This may be because they are waiting to accumulate sufficient capital to make further 
stock purchases; or they tend to monitor their portfolios infrequently; or they are simply scared 
of choosing the wrong investments; or they feel that as long-term investors, they can defer action; 
or they put off the time-consuming research necessary to make informed portfolio choices. Long 
uninterrupted periods of inactivity may be viewed as a form of investor inertia. 

2.4.1 Evidence of inertia 

Investor inertia is a common experience and is well documented. The New York Stock Exchange 
(NYSE)'s survey of individual shareownership in the United States, "Shareownership2000" |33j, 
demonstrates that many investors have very low levels of trading activity. For example they find 
that "23 percent of stockholders with brokerage accounts report no trading at all, while 35 percent 
report trading only once or twice in the last year" . The NYSE survey also reports (Table 28) that 
the average holding period for stocks is long, for example 2.9 years in the early 90's. Empirical 
evidence of inertia also appears in the economic literature. For example, Madrian and Shea |67j 
looked at the reallocation of assets in employees' individual 401 (k) (retirement) plans and found 
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"a status quo bias resulting from employee procrastination in making or implementing an optimal 
savings decision." A related study by Hewitt Associates (a management consulting firm) found that 
in 2001, four out of five plan participants did not do any trading in their 401(k)s. Madrian and 
Shea explain that "if the cost of gathering and evaluating the information needed to make a 401 (k) 
savings decision exceeds the short-run benefit from doing so, individuals will procrastinate." The 
prediction of Prospect Theory (see [56]) that investors tend to hold onto losing stocks too long has 
also been observed in [77]. Another typical cause is that small investors seem to find it difficult to 
reverse investment decisions, as is discussed even in the popular press. A recent newspaper column 
(by Russ Wiles in the Arizona Republic, November 30, 2003) states: "Perhaps more than anything, 
investor inertia is a key force (in financial markets). When the news turns sour, people tend to 
hold off on buying rather than bail out. In 2002, the toughest market climate in a generation and 
a year with ample Wall Street scandals, equity funds suffered cash outflows of just one percent." 

2.4.2 Inertia and long range dependencies in financial time series 

One of the outcomes of a limit analysis of an agent-based model of investor inertia is a stock 
price process based on fractional Brownian motion, which exhibits long-range dependence (that 
is correlation or memory in returns). This is discussed in Section 13.11 In particular, the limit 
fluctuation process is a fractional Brownian motion. 

We recall that fractional Brownian motion B with Hurst parameter H G (0, 1] is an almost 
surely continuous and centered Gaussian process with auto-correlation 

E {B*B?} = \ (\t\ 2H + \s\ 2H -\t- s\ 2H ) . (2.6) 

Remark 2.1 Note that the case H = ^ gives standard Brownian motion. Also note that the 
auto- correlation function is positive definite if and only if H £ (0, 1]. 

Bayraktar et al. [8] studied an asymptotically efficient wavelet-based estimator for the Hurst 
parameter, and analyzed high frequency S&P 500 index data over the span of 11.5 years (1989- 
2000). It was observed that, although the Hurst parameter was significantly higher than the efficient 
markets value of H = ^ up through the mid-1990s, it started to fall to that level over the period 
1997-2000 (see Figure[T]). This might be explained by the increase in Internet trading in that period, 
which is documented, for example, in NYSE's "Shareownership2000" [33] . Barber and Odean [2], 
and Choi et al. [22], in which it is demonstrated that "after 18 months of access, the Web effect 
is very large: trading frequency doubles." Indeed, as reported in [3], "after going online, investors 
trade more actively, more speculatively and less profitably than before" . Similar empirical findings 
to that of [8] were recently reached, using a completely different statistical technique by Bianchi 

i- 

Thus, the dramatic fall in the estimated Hurst parameter in the late 1990s can be thought of 
as a posteriori validation of the link the limit theorem in [5] provides between investor inertia and 
long-range dependence in stock prices. We review this model in Section 13.11 An extension based 
on state dependent queuing networks with semi-Markov switching is discussed in Section 13.21 
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Figure 1: Estimates of the Hurst exponent of the S&P 500 index over 1990s, taken from Bayraktar, Poor 
and Sircar JSjj. 

3 Microstructure Models with Inert Investors 

We illustrate the use of microstrucure, or agent-based models, combined with limit theorems by 
focusing on investor inertia as a very common characteristic among small and casual market partic- 
ipants. In Section [3.11 we summarize earlier work [5] that established a mathematical link between 
inertia, long-range dependence in stock returns and potential short-lived arbitrage opportunities 
for other 'sophisticated' parties. Section 13.21 contains an extension allowing for feedback effects 
from current prices into the agents' order rates. 

3.1 A Microstructure Model without Feedback 

We now introduce the basic concepts and notation of the market microstructure model analyzed 
in [5] that will serve as basis for the more sophisticated model in Section I3.2L We start with a 
financial market with a set A := {a±, «2, • • • > Qjv} of agents trading a single risky asset. Each agent 
a G A is associated with a continuous-time stochastic process x a = {xf }t>o on a finite state space 
E describing his trading activity. 

We take a pragmatic approach to specify the demand. Instead of formulating an individual 
optimization problem under budget constraints for the agents, we start right away with the agent's 
order rates. The agent a £ A accumulates the asset at a rate ^jx" at time t > 0. Here xf may 
be negative indicating that the agent is selling. The random process VP = {^t}t>o describes the 
evolution of the size of a typical trade. It can also be interpreted as a stochastic elasticity coefficient 
(the reaction of the price to the market imbalance) . We assume that ^ is a continuous non- negative 
semi-martingale which is independent of the processes x a and that S E. The agents do not trade 
at times when x% = 0. The holdings of the agent a G A and the "market imbalance" at time t > 
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are thus given by, respectively, 

I ^ s x a s ds and V f ^ s x a s ds. (3.1) 

Remark 3.1 In our continuous time model, buyers and sellers arrive at different points in time. 
Hence the economic paradigm that a Walrasian auctioneer can set prices such that the markets clear 
at the end of each trading period does not apply. Rather, temporary imbalances between demand 
and supply will occur. Prices will reflect the extent of market imbalance. 

All the orders are received by a single market maker. The market maker clears all trades 
and prices in reaction to the evolution of market imbalances, the only component driving asset 
prices. Reflecting the idea that an individual agent has diminishing impact on market dynamics 
if the number of traders is large, we assume that the impact of an individual order is inversely 
proportional to the number of possible traders: a buying (selling) order increases (decreases) the 
logarithmic stock price by 1/N. The pricing rule for the evolution of the logarithmic stock price 
process S N = {S^}t>o is linear and taken to be: 

dS t = ^J2^dt. (3.2) 

aeA 

In order to incorporate the idea of market inertia, the agents' trading activity is modelled by 
independent and identically distributed semi-Markov processes x a . Semi-Markov processes are 
tailor-made to model individual traders' inertia as they generalize Markov processes by removing 
the requirement of exponentially distributed, and therefore thin-tailed, holding (or sojourn) times. 
Since the processes x a are independent and identically distributed, it is enough to specify the 
dynamics of some "representative" process x = {x} t >o- 

3.1.1 Semi- Markov Processes 

A semi-Markov process x defined on a probability space P) is specified in terms of random 

variables £ n : Q — > E and T n : Q, — ► R + , satisfying = T\ < T\ < ■ ■ ■ almost surely and 

P{£n+i = j,T n+ \ -T n < t|£i, ...,£ n ;Ti, ...,T n } = P{£ n+ i = j,T n+ i - T n < t\^ n } 

for each n G N, j G E and all t G E+, through the relation 

x t = J2^ nl lTn,T n+1 )(t)- (3.3) 

n>0 

In economic terms, the representative agent's mood in the random time interval [T n ,T n+ i) is given 
by £ n . The distribution of the length of the interval T n+ \ — T n may depend on the sequence {£ n } n gN 
through the states £ n and £ n +i. This allows us to assume different distributions for the lengths of 
the agents' active and inactive periods, and in particular to model inertia as a heavy-tailed sojourn 
time in the zero state. 
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Remark 3.2 In the present analysis of investor inertia, we do not allow for feedback effects of 
prices into agents ' investment decisions. While such an assumption might be justified for small, 
non-professional investors, it is clearly desirable to allow active traders' investment decisions to be 
influenced by asset prices. We discuss such an extension in the next section. 



We assume that x is temporally homogeneous under the measure P, that is, 
Q(i,3, t) = P{£n+i = 3, T n+1 -T n < t\£ n = i} 



(3.4) 



is independent of n G N. By Proposition 1.6 in [17] . this implies that {£ n }neN is a homogeneous 
Markov chain on E whose transition probability matrix (pij) is given by 

Pij = lim Q(i,j,t). 

t— >oo 

Clearly, x is an ordinary temporally homogeneous Markov process if Q takes the form 

Q(i,j,t)= Pij (l-e~ Xit y (3.5) 

We also assume that the embedded Markov chain {£ n }neN satisfies P ij > so that {£ n }neN has 
a unique stationary distribution. The conditional distribution function of the length of the n- 
th sojourn time, T n+ \ — T n , given £ n +i and £ n is specified in terms of the semi-Markov kernel 
{Q(i, j,t);i, j EE, t > 0} and the transition matrix P by 



G(i,j,t) := Q( - h3 > l) = P{r n+1 - T n < t\£ n = i, Cn+i=j}. 
Pij 

The semi-Markov processes are assumed to satisfy the following conditions. 
Assumption 3.3 (i) The average sojourn time at state i 6 E is finite: 

mi := E[T n+1 - T n |£ n = i] < oo. 
Here E denotes the expectation operator with respect to P. 

There exists a constant 1 < a < 2 and a locally bounded function L : I 
slowly varying at infinity (e.g. log), i.e., 

L{xt) 



(3.6) 



(3.7) 



which is 



lim 



1 



for all 



x > 0, 



t^oo L(t) 

such that 

P{T n+1 -T n > t\U = 0} ~ t- a L(t) 

Here we use to notation f(t) ~ g(t) for two functions f,g 
lim^oo f{t)/g{t) = 1. 

(Hi) The distributions of the sojourn times at state i ^ satisfy 

P{Tn+l-T n >t\£ n = i} 



(3.8) 

to mean that 



lim 



t-(«+i)L(t) 



0. 



(3.9) 
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(iv) The distribution of the sojourn times in the various states have continuous and bounded 
densities with respect to Lebesgue measure on M + . 

The key parameter is the tail index a of the sojourn time distribution of the inactive state 
zero. Condition (|3.8p is satisfied if, for instance, the length of the sojourn time at state G E 
is distributed according to a Pareto distribution. The idea of inertia is then reflected by (|3,9p : 
the probability of long uninterrupted trading periods is small compared to the probability of an 
individual agent being inactive for a long time. In fact, it is natural to think of the sojourn times in 
the various active states as being thin tailed as in the exponential distribution since small investors 
typically do not trade persistently. 

3.1.2 A Limit Theorem for Financial Markets with Inert Investors 

We assume that the semi-Markov processes x a are stationary. Stationarity can be achieved by a 
suitable specification of the common distribution of the initial states and initial sojourn times. We 
denote the resulting measure on the canonical path space by P*. Independence and stationarity 
of the semi-Markov processes guarantees that the logarithmic price process can be approximated 
pathwise by the process {s t }t>o defined by 



St = n ^> s ds where \x := E*2q 
Jo 



10 

when the number of agents grows to infinity. Our functional central limit theorem for stationary 
semi-Markov processes shows that after suitable scaling, the fluctuations around (st)t>o can be 
approximated in law by a process with long range dependence. The convergence concept we shall 
use is weak convergence with respect to the measure P* of the Skorohod space B of all right 
continuous processes. We write C-liuin^oo Y n = Y if {Y n } ng N is a sequence of D-valued stochastic 
processes that converges weakly to the process Y. 

The convergence result is formulated in terms of a scaling limit for the processes {x^ t }t>o 
(TEN). For T large, "speeded-up" semi-Markov process. In other words, the investors' 

individual trading dispensations are evolving on a faster scale than \P. Observe, however, that we 
are not altering the main qualitative feature of the model: agents still remain in the inactive state 
for relatively much longer times than in an active state. In the rescaled model the logarithmic asset 
price process S N,T is given by 



i I ^uX a Tu du. (3.10) 

J ° aeA 



The central limit theorem allows us to approximate the fluctuations around the first order approx- 
imation as N — > oo. In terms of the Gaussian processes X T and Y T defined by 

i N r l 

If^-lim T 1 -^^-^) and Y t T 4 Xjds, (3.11) 

N->co V N ~i JO 

with H = (3 — ct)/2, the fluctuations around the first order approximation can be approximated 
by an integral of the elasticity coefficient with respect to Y T : 

r-t 



C- lim Vn{s?' T -fit} ={ I ^sdY^X 

N~+oo I )0<t<l [Jq J 



0<t<l 
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In order to see more clearly the effects of investor inertia, we rescale the price process in space 
and time and T tend to infinity. In a benchmark model with many agents where VP = 1 these, 
fluctuations when suitably normalized, can be described by a fractional Brownian motion B H if 
T — > oo. The Hurst coefficient is related to the degree of investor inertia. 

Theorem 3.4 ([5]) Let H = ^jp. Assume that *S> = \, that Assumption ] 3. 3\ holds and that /i / 0. 
Then there exists a > such that 



C- lim C- lim r i-g l^_ \s?' T -nt\ = {aB? k (3.12) 

To generalize this result to a market in which the agents' order rates are coupled by a stochastic 
elasticity coefficient as in (|3.2p . we need the following approximation result for stochastic integrals 
of continuous semi-martingales with respect to fractional Brownian motion. 

Theorem 3.5 (J^j) Let {* n }neN be a sequence of good semimartingales and {Z n } n ^ be a sequence 
ofH)-valued stochastic processes that satisfy 

(i) The sample paths of the processes Z n are almost surely of zero quadratic variation on compact 
sets, andF{Z$ = 0} = 1. 

(ii) The stochastic integrals J fy n dZ n and J Z n dZ n exist as limits in probability of Stieltjes-sums, 
and the sample paths t J Q ' Z™dZ2 and t \— > L ^™dZ2 are cadlag. 

If ^ is a continuous semimartingale and if B H is a fractional Brownian motion process with Hurst 
parameter H > ^> then the convergence £-lim n ^ 00 (^I' ra , Z n ) = (^f,B H ) implies the convergence 

C- lim f ^ n ,Z n , [ $> n dZ n ) = (v,B H , [ ^dB H 



As an immediate corollary to Theorem 13.51 we see that the fluctuations of the price process 13.101 
around its first order approximation converge in distribution to a stochastic integral with respect 
to fractional Brownian motion. 

Corollary 3.6 Let ^ be a continuous semi martingale with Doob-Meyer decomposition ^ = M+A. 
If E{[M, M]^} < oo, E{|^4|^} < oo and n^O, then there exists a > such that 

C- lim C- lim T 1 - H ^==is?' T - n C ^ s ds\ ={a[\ s dB^\ . (3.13) 
r-oo W-oc s/W)\ Jo ) <t<i I Jo )o<t<i 

The increments of a fractional Brownian motion with Hurst coefficient H € 1] are positively 
correlated. The correlation increases in H. Thus, the limit theorem reveals that, in isolation, 
investor inertia may lead to long range dependence in asset returns. Indeed, a greater degree 
of inactivity, represented by a smaller tail index a, leads to a larger H, and so greater positive 
correlation between returns. Since fractional Brownian motion is not a semimartingale, it may 
also lead to arbitrage opportunities for other traders whose impact has not been considered in the 
model so far. Explicit arbitrage strategies for various models were constructed in, e.g. [6]. 
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Remark 3.7 In a model without inertia where all the sojourn time distributions are thin-tailed, 
the logarithmic stock price fluctuations can be approximated in law by a process of the form 

\lv s dW s \ (3.14) 
Uo J o<t<i 

where W is a standard Brownian motion. Thus, when all traders' mood processes are standard 
Markov processes and \E' is constant, we recover in the limit the standard Black- Scholes-Samuelson 
geometric Brownian motion model. 

The approach of studying queuing systems through their limiting behaviour has a long history 
in many applications, see [83J, for example. This analysis of investor inertia built upon the works of 
Taqqu et al. [80\ on internet traffic. However, even the simple model we have discussed so far shows 
how economic applications lead to new mathematical challenges: in the teletraffic application, it is 
sufficient to consider a binary (on/off) state space, but when agents buy, sell or do nothing, there 
must be at least three states. This requires different techniques from the binary case. Our functional 
central limit theorems for stationary semi-Markov processes may also serve as a mathematical basis 
for proving heavy-traffic limits in the multilevel network models studied in, e.g. [30J and [29|. 

3.2 A Limit Theorem with Feedback Effects 

The model in the previous section assumes that investors' actions affect the price, but prices did 
not affect the agents' demands. This assumption might be justified for Internet or new economy 
stocks where no accurate information about the actual underlying fundamental value is available. 
In such a situation, price is not always a good indicator of value and is often ignored by unin- 
formed small investors. In general, however, it is certainly desirable to allow for feedback effects 
from current prices into the agents' order rates. In this section we extend our previous model to 
allow for feedback effects from prices into the agents' order rates. At the same time we provide 
a unified mathematical framework for analyzing microstructure models with asynchronous order 
arrivals. Our approach is based on methods and techniques from state dependent Markovian ser- 
vice networks. Mathematically, it extends earlier results in [1] beyond semi-Markov models with 
thin-tailed sojourn time distributions. 

3.2.1 The dynamics of logarithmic asset prices 

Let us now be more precise about the probabilistic structure our model. We assume that the agents' 
orders arrive with an order rate that depends on the price and the investor sentiment. Each order 
is good for one unit of the stock. Specifically, we associate to each agent a £ A two independent 
standard Poisson processes {II" (*)} t>0 an d {rP (*)} t > > a stationary semi-Markov process x a on 
E satisfying Assumption 13.31 and bounded Lipschitz continuous rate functions \± : E x R — ► R + . 
The rate functions along with the Poisson processes ILj_ specify the arrivals times of buying and 
selling orders. The agent's holdings at time t > are given by 

II« (J* A + «, S») dv\ - U a _ A_ «, S») dv\ (3.15) 
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where {S^ }t>o denotes the logarithmic asset price process. As before, a buying (selling) order 



increases (decreases) the logarithmic price by 1/N. Assuming for simplicity that Sq = 0, we thus 



obtain 

s ? = Tr E n + ( f A + s ") du )-JjH n - ( f x - «> S u) du) . (3.16) 

Remark 3.8 (i) In the model studied in the previous section, the agents continuously accumu- 
lated the stock at rates specified by semi-Markov processes. Our current models assume that 
stocks are purchased at random points in times. The arrival times of buying and selling times 
follow exponential distributions conditional on random arrival rates that depend on current 
prices and exogenous semi-Markov processes. 

(i) As before, we think of x a as being the investor's "mood" (for trading) process. Loosely speak- 
ing, X + (xf, s)—X-(x1, s) may be viewed as the agent's excess demand at time t at a logarithmic 
price level s, given his trading mood x%. 

(ii) To develop a model of interaction, in which the participants are inert, out of 113. 15\) . it is 
natural to assume that A±(0, s) = and that the buying and selling rates A+(x, •) and A_(x, •) 
are increasing, rep. decreasing, in the second variable meaning that meaning high (low) prices 
temper buying (selling) rates. 

The sum of independent Poisson processes is a Poisson process with intensity given by the sum 
of the intensities. As a result, the logarithmic price process satisfies the equality 

\7 n + (E jf A + «> S u) " (E jf A - «> S u) ^ ( 3 - 17 ) 

in distribution where n + and II_ are independent standard Poisson processes. Since our focus 
will be on a limit result for the distribution of the price process as the number of agents grows to 
infinity, we may with no loss of generality assume that the logarithmic price process is defined by 
(pTTTft rather than (pTTBj) . 

Assumption 3.9 1. The rate functions X± are uniformly bounded. 

2. For each x £ E, the rate functions X±(x, •) are continuously differentiate with first derivative 
bounded in absolute value by some constant L. 

Our convergence results will be based on the following strong approximation result which allows 
for a pathwise approximation of a Poisson process by a standard Brownian motion living on the 
same probability space. 

Lemma 3.10 (1 6 If ) A standard Poisson process {H(t)} t>0 can be realized on the same probability 
space as a standard Brownian motion {B(t)}t>o in such a way that the almost surely finite random 
variable 

\U(t) - t - B(t)\ 

SUP : ; 

t>o log(2 V t) 

has a finite moment generating function in the neighborhood of the origin and in particular finite 
mean. 
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In view of Assumption 13.91 (i) , the strong approximation result yields the following alternative 
representation of the logarithmic asset price process: 



S t = Jf IE jf A & S ") du + B + (E jf A + & S u) A 



(3.18) 

logiV N 



iV 

where A(x", •) denotes the excess order rate of the agent a £ A, given his mood for trading 
and O (log N/N) holds uniformly over compact time intervals. Using this representation of the 
logarithmic price process our goal is to prove approximation results for the process {S^}t>o- In 
a first step we show that it can almost surely be approximated by the trajectory of an ordinary 
differential equation ("fluid limit"). In subsequent step, we apply a result from [5] to show that, 
after suitable scaling, the fluctuations around this first order approximation can be described in 
terms of a fractional process {Z t }t>o of the form 

dZ t = fitZtdt + (JtdBf 1 . 

In a benchmark model without feedback, where the order rates do not depend on current prices, 
the process {Z t }t>o reduces to a fractional Brownian motion. That is, we recover the type of results 
of Section 13.1.21 with the alternative model presented in this section. 

3.2.2 First order approximation 

In order to prove our first convergence result, it is convenient to denote by 

X(x, s) = X+(x, s) — A_(x, s) (3.19) 

the accumulated net order rate at a given logarithmic price level s£l and trading mood x G E 
and by 

A( S ) = A + ( S )-A-( S ) 

the expected excess order flow where 

A±(s) = / \±(x, s)v(dx), 
Je 

and v is the stationary distribution of the semi-Markov process xt- We are first going to show 
that in a financial market with many agents the dynamics of the logarithmic price process can be 
approximated by the solution {st}t>o to the ODE 

j t s t = \(s t ), (3.20) 

with initial condition sq = 0. To this end, we need to prove that the average excess order rate 
converges almost surely to the expected excess order flow uniformly on compact time intervals. 
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Lemma 3.11 Uniformly on compact time intervals 

i N 
im — 

N 



I N ft ft _ 

lim — :) / A±«,s u )du = / \±(s u )du F*-a.s. (3.21) 



Proof: The stationary semi Markov processes x a are independent, and so the random variables 
Jo s u )du (a = 1, 2, ...) are also independent. Thus, the law of large numbers for independent 
random variables along with Fubini's theorem (to exchange the sum and the integral) and bounded 
convergence theorem (to exchange the limit and the integral) yields convergence for each t. In order 
to prove that the convergence holds uniformly over compact time intervals we will use uniform law 
of large numbers of [71]. Denoting D^[0,t] the class of all cadlag functions y : [0,i] — > E we need 
to show that the maps q± : De[0, t] x [0, t] — ► K defined by 

Q±(y,t) = X±(y(u),s u )du 
Jo 

are continuous. Since the rate functions are bounded, it is enough to show that the map y i— > 
Jo X±(y(u), s u )du is continuous uniformly over compact time intervals. 

To this end, we denote by d the metric defined in (3.5.2) in [31] which induces the Skorohod 
topology in De[Q, t] and recall that hnin^oo d(y n ,y) = if and only if 

lim sup \y n o r n (s) - y(s)\ = (3.22) 

n ^°°0<s<t 

for a suitable sequence of strictly increasing time-shifts r n ; see [34] page 117 for details. Let {y n } 
denote a sequence in De[0, t] that converges to y and put 

\±{u) = X±(y n (u),s u ). 

In view of the transformation formula for Lebesgue integrals and because r(0) = and t" 1 ^) < t 
we obtain 

ft f-n 1 ^) ft 

/ [\±(u) - \±(u)] du = [XI o T n {u)T' n {u) - X±(u)] du - / X ± {u)du 

Jo Jo Jt" 1 ^) 



TnHt) 

[A J o T n (u) — A-t(u)] du 



r T n 1 (t) ft 
+ / A^ o T n {u)[T' n (u) - l]du - / X±(u)du. 

JO Jr-Ut) 



By (3.5.5)-(3.5.7) in 



lim sup \T n {u) — 1 1 = and lim sup \r n (u) — u\ = 

n ^°°o<M<t n ^°°o<M<t 

so that the last two terms on the right hand side of the inequality above vanish uniformly on 
compact time intervals. As far as the first term is concerned, observe that boundedness of the rate 
function's derivative with respect to the second argument yields 

\X±(y n ° r n (n),s Tn(n )) - X±(y(u),s u )\ < L \y n o r n («) - y(u)\ + L\s o r n (u) - s(u)\. 
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As a continuous function s is uniformly continuous over compact time intervals. This, along with 
(pT22j) yields 

lim sup \X±(y n oT n (u),s T f u) )-\±(y(u),s u )\=Q 
n -^°°o<«<t 

so that the maps q± are indeed continuous. Thus, the uniform law of large numbers yields 

N 



1 N 
lim — > 



0=1 



1 -> l' u 

lim "1^7 E / \±(Xv,s v )dv = \±(n,s u ) 

N— >oo IV 1 — * /n 

a=l Jy> 



□ 



almost surely uniformly on compact time intervals. 

We are now ready to state and prove our functional law of large numbers. 

Theorem 3.12 As N — > oo ; the sequence of stochastic processes {S^}t>o (N € N) converges 
almost surely to the deterministic process {st}t>o-' 



lim 5+ = st 

N->oc 1 



-a.s. 



where the convergence is uniform over compact time intervals. 

PROOF: In view of the strong approximation result formulated in Lemma 13.101 and because the 
rate functions are uniformly bounded, 



n± E/ A ± {<,Su)du) ~E / *±{<,S*[)du-B± E/ A ± 



is of the order O (logiV) almost surely where B± are the Brownian motions used in (|3,18p . Since 
the rate functions are uniformly bounded, the law of iterated logarithm for Brownian motion yields 



lim snpli? ± f E r A± {x a v ,S») dv) = 
n-*oo u < t J\ Jo J 



-a.s. 



It follows from this and Lemma 13.111 above, that the quantities 



B + (E jf A + «> dn ) " B - (E / A - «> S u) du 



and 



Af 4 



N r t 



I f 

atE / - A(s u )}dn 

o=l ^° 



converge to zero uniformly over compact time intervals as N — > oo. 

Let us now fix e > 0. Due to Lemma 13.101 there exists iV* € N such that for all N > N* and 
uniformly on compact time sets, for I < t we can write 

N rl 

i)du 



< 



^E / A «,^)du- / A(s„ 
7V «=i Jo 

-E / {\(xlS%)-\(x a u ,s u )}du 
7V «=i ^ 



-a.s. 
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Lipschitz continuity of the rate functions yields 



ISf-sA < l[ sup \S?-s r \du + Af + Bf f + e 

JO 0<r<u 

< L f sup \S?-s r \du+ sup sup B? + e P*-a.s. 

JO 0<r<u 0<r<t 0<r<t 



for some L > and so 



sup |<S^ — s r \ <L sup |<S^ — s r \du+ sup + sup B^ + e P*-a.s. (3.23) 

0<r<t Jo 0<r<« 0<r<t 0<r<t 

Now, an application of Gronwall's lemma yields 

sup |Sf - s r \ < ( sup + sup B? + e J e L * P*-a.s. 

0<r<i V0<r<t 0<r<t / 

for all N > N* . This proves our assertion. □ 



3.2.3 Second order approximation 

In this section we analyze the fluctuations of the logarithmic price process around its first order 
approximation. We are interested in the distribution of asset prices around their first order ap- 
proximation as N —* oo. In view of the representation (|3.18p and by self-similarity of Brownian 
motion we may thus assume that {S^}t>o is defined by the integral equation: 




(3.24) 



As we shall see, the fluctuations around the first order approximation are driven by two Gaussian 
processes. The first, 

X t = B + (j % \+{s u )dv\ - B^ (J* \„(s u )dv\ , (3.25) 

captures the randomness in the agents' trading times. The second, {Yt}t>o, is defined in terms of 
the integral of a non-stationary Gaussian process whose covariance function depends on the first 
order approximation. It captures the second source randomness generated by the agents' trading 
activity. Specifically, 

Y t = f y s ds, (3.26) 
Jo 

where {yt}t>o denotes the centered Gaussian process whose covariance function 7 is given by the 
covariance function of the stochastic process {X(xt,st)}t>o, i-e., 

7(t, u) = E[A(xt, s t )X(x u , s u )] - X(s t )X(s u ). (3.27) 
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It turns out that the fluctuations can be approximated in distribution by the process {Zt}t>o which 
satisfies the integral equation 

Z t = f \'{s u )Z u du + Y t +X t . (3.28) 
J o 

Our goal is to establish the following second order approximation for the asset price process in 
an economy with many market participants. 

Theorem 3.13 The fluctuations of the market imbalance {S^}o<t<i around its first order approx- 
imation can be described by the process {^}o<t<i defined in \3.28\) . More precisely, 

C- lim ^{Sf-stj^^iZt}^. 



The proof of Theorem 13.131 requires some preparation. For notational convenience we introduce 
stochastic processes Q N = {Q^oKt^i, Y N = {Y t N }o< t <i and X N = {Xj^}o< t <i by, respectively, 

Qf A ^n {S ? - s t ) and Y t N 4 V f X (x a u ,s u )- X(s u ) ^ ^ 

«=i Jo v N 

and 

X? = B + (± £ jf A + «, S?) duj - B_ f 1 g jT A_ «, 5^ ^ . (3.30) 



We first prove convergence in distribution of the sequence {(X , Y )}jveN to (X,Y). 

Proposition 3.14 The sequence {(X , Y n )}n^ converges in distribution to the process (X, Y) 
defined by $3.25\) and (EI 



Proof: For any a € (0, |) and T > 0, there exist integrable and hence almost surely finite random 
variables M± such that for all t\, t% < T we have 

\B±(ti) - B±(t 2 )\ < M±\h - t 2 \ a P*-a.s., 

see, for instance, Remark 2.12 in [57]. Thus, the first order approximation shows that the sequence 
of processes {X^jjvgN converges almost surely to X on any compact time interval. Since the 
processes 

r A(a;°, s u ) — X(s u ) 
= du 

have Lipschitz continuous sample paths and the semi-Markov processes are independent, the central 
limit theorem for Lipschitz processes Q83J, Corollary 7.2.1) shows that {Y^jveN 

converges m 

distribution to the Gaussian process Y. As a result, both sequences {X^j-TveN and {F^ljveN are 
tight. Since {X^jngM is also C-tight, the sequence {(X , Y^jiveN is tight. It is therefore enough 
to prove weak convergence of the finite dimensional distributions of the process (X N ,Y N ) to the 
finite dimensional distributions of (X, Y). 
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In order to establish weak convergence of the one-dimensional distributions we fix a Lipschitz 
continuous functions with compact support F : M. 2 — > M. We may with no loss of generality assume 
that both the Lipschitz constant and the diameter of the support of F equal one. In this case 

J F (X t N , Y t N ) dF* - J F (X t , Y t N ) dF* < J mm{\X t N - X t \,l}dF*. 

In view of the convergence properties of the sequence {X^jjveN) there exists, for any e > 0, a 
constant N* G N such that 

sup / mm{|Zf - X t \,l}dF* < e for all N > N* . 

0<Kli 



This yields 



lim 

N— >oo 



J F (X t N ,Yf) dF* - J F(X t ,Y t N ) 



dF* 



Since the random variables Xt and Y t are independent, we also have that 



lirn^ J F(X t , Y t N ) dF* = J F {X u Y t ) dF* . 

This proves vague convergencj§ of the one-dimensional marginal distributions of (X N ,Y N ) to the 
one-dimensional distributions of (X, Y) and hence weak convergence. Weak Convergence of the 
finite dimensional distributions follows from similar considerations. □ 

The following "compact containment condition" is key to the second order approximation. 

Lemma 3.15 (i) The sequence of stochastic processes {Q n }n&n ^ bounded in probability. That 
is, for any e > 0, there exists N* € N and K < oo such that 



sup | | > K 

0<t<l 



< e for all N > N* . 



(3.31) 



(ii) If f = {/t }t>o be a sequence of non-negative random processes such that 

»l 



then, for all 5 > 0, 



lim 

N^>oo 



lim P* 

N— >oo 



fudu = in probability, 



sup 

o<a<i 



Qufudu 



> 5 



0. 



(3.32) 



Proof: 



A sequence of probability measure {/j„} converges to a measure [J, in the vague topology if lim, woo f fdfi n = f fd^i 
for all continuous functions / with bounded support. The vague limit /j, is not necessarily a probability measure. 
However, if there is an a priori reason that n is a probability measure, then weak convergence of {/J. n } to \x can be 
established by analyzing integrals of continuous and hence Lipschitz continuous functions with bounded support. See 
e.g. g] and [TO]. 
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(i) The strong approximation for Brownian motion yields the representation 



Q 



JtEa=l{^«,S^)-\(x a u , Su )}du 



+ Y t N + X t N + O 



logiV 



(3.33) 



By Proposition 13.141 the sequence {(X N ,Y N )} n& ^ is tight, and hence it is bounded in prob- 
ability (see e.g. [29]). As a result, Lipschitz continuity of the rate functions yields 

sup \Qf\<L [ sup \Qu\du+ sup \Y t N \ + sup \Xf | + O ( lo " A 

0<t<l JO 0<t<u 0<t<l 0<t<l 

for some L > 0. Hence, by Gronwall's inequality, 



sup |Qf | < e iT 

0<<<1 



sup \Y t N \ + sup \X t N \ + o(^j£ 
o<t<i o<t<i \ vjV 



-AT I 



-a.s. 



This proves (i). 

(ii) Let us fix e > 0. There exists a constant iV* such that when ./V > N* there exist sets and 
An such that 

f 1 e e 

/ fudu<- on^Tv and such that F*[Q N ] > 1 - -. 
Jo 2 2 



and 



Hence 



sup | | < K on A N and such that P*L4at] > 1 - -. 
o<t<i 2 



sup 

0<t<l 



I Qufudu < sup / f*du<Ke onA N nn N . 

Jo o<t<i Jo 



□ 



Proof of Theorem \3.13[ Let us first define a sequence of stochastic processes Q = {Qt}o<t<i by 

Jo 

By the continuous mapping theorem and Lemma 13.141 the sequence {Q^jjveN converges in distri- 
bution to the process Z defined in (|3.28p . It is now enough to show that 



lim sup \Q t — Q t | = in probability. 

iV^oo <t<l 



(3.34) 



To this end, let Ef = — . ^From the definition of and the representation (|3,33p of Q t 
we obtain 



N 



\'(s u )E N du+^= [ J2\ X K S u) ~ \«,s u ))du 
VN Jo a=1 I J 

jf \'( Su )E^du + J*(±J^ A '«> s «) " A '( S «)J Qu du 



X(s u )Q u du 



25 



The second equality follows from the mean value theorem for A(x", •), 



N 



where ^ lies between jtS^ and s u . We put 



N N 



a=l 



a=l 



in order to obtain 



sup \E?\<L / sup \E?\du + 

0<s<t JO 0<s<u 



0<s 



+ 


sup / 




0<s<t Jo 



The processes |/ ' | and |/ ' | satisfy the condition (I3.32p of Lemma 13.151 by the law of large 
numbers. Thus, an application of Gronwall's lemma yields (|3.34p . □ 

3.2.4 Approximation by a fractional Ornstein-Uhlenbeck process 

So far we have shown that the fluctuations of the logarithmic price process around its first order ap- 
proximation can be described in terms of an Ornstein-Uhlenbeck process Z driven by two Gaussian 
processes X and Y. In order to see more clearly the effects of investor inertia on asset processes 
we need to better understand the dynamics of Y. As before, this will be achieved by a proper 
scaling of of the semi-Markov processes x a in time and the price process in space. Specifically, we 
introduce a family of processes S N ' T (T G N) with initial value by 



S, 



N.T 



The strong approximation result for Poisson processes with respect to Brownian motion allow us 
to represent the process {Sf ,T } t >o as in (|3.18p with the semi-Markov processes {x® } t >o replaced 
by the "speeded-up" processes {x^ t } t >o- Moreover, by Lemma f3. Ill the sequence of processes 



A 



N,T A 



1 



N 



a=l JU 



{\(xt u ,s u ) - X(s u )}du 



converges to zero uniformly over compact time intervals as N —* oo. Following the same line of 
arguments as in the proof of Proposition 13.121 it can then be shown that for any T > 



lim b f ' = St 

N->oo 



-a.s. 



(3.35) 



Here {st}t>o denotes the deterministic process defined by the ordinary differential equation (|3.20p 
with initial condition sq = and the convergence holds uniformly over compact time intervals. 
Thus, the first order approximation is independent of T. By analogy to (|3.24|) - (|3.28j) introduce a 
Gaussian process Y T by 



T A 



Vs ds ' 



(3.36) 
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where {yf}t>o denotes the centered Gaussian process with covariance function 

-y T (t, it) = E[X(xTt, s t )X(x T u, s u )] - X(s t )X(s u ). 

Following the same arguments in the proof of Theorem 13,131 we see that as the number of agents 
tends to infinity the price fluctuations round the fluid limit can be approximated in distribution by 
a process {Zj} t >o of the form 



Zj = J ~X'{s u )Zldu + Y t T + -±=X t . 



Proposition 3.16 For any T, the fluctuations of the logarithmic price process {S t ' }o<t<i around 
its first order approximation can be described by the process {Z^}o<t<i- More precisely, 

C- lim VN\S?' T - s t \ = {Z?} Q<t<1 . 

N^oo I 1 J 0<i<l L ' 

To take the T-limit, we need the following assumption on the structure of the rate functions. 
Assumption 3.17 The rate function X defined in \3. 19\) can be written as 

X(x,s) = f(x)g(s) + h(s). (3.37) 
Moreover, the function f in is one-to-one and fi = /(0) ^ E*/(xo). 

Example 3.18 The previous assumption is always satisfied if (xf)t>o is a stationary on/off pro- 
cess, i.e., if E = {0,1}. In this case 

a _ X(x<Z,s t )-X(0,s t ) 
1 A(l, St )-A(0, St ) ' 



and the representation (3.31) holds with 

f(x) = x, g(s) = A(l, s) — A(0, s) and h(s)=X(0,s). 

We are now ready to show that the fluctuations of the logarithmic stock price around its first 
order approximation behaves like a fractional Ornstein-Uhlenbeck process. 

Theorem 3.19 Under the Assumptions \3.9\ and \3.17\ we have that 



C- lim C- lim T l - H 1 {s?' T - s t \ = {Z t } 0<t<1 

Here Z denotes unique solution to the stochastic differential equation 

dZ t = X'(s t )Z t dt + ag{s t )dB? 

where B is a fractional Brownian motion with Hurst coefficient H = The integral with 

respect to B H is understood as a limit in probability of Stieltjes sums. 
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Proof: The proof uses modifications of arguments given in the proof of Theorem 13.131 and the 
approximation result for integrals with respect to fractional Brownian motion in [5]. 

(i) In a first step we study the dynamics of the process {Y t ' }(>o defined by 

y N,T = /"* K^Tu^u) - A(g u ) du 
1 L . /at 

a=l 

Under Assumption 13. 171 we can write 



o VN 



Y t N,T = E / "7X7 WtuWsu) + h{s u ) - X(s u )] du 
a=1 Jo VN 

= Yl / ~7tt - A] g(s u )du 

a=1 Jo VN 



(3.38) 



Since / is one-to-one, (f{x^))t>o is a semi-Markov process that has the same sojourn time 
structure as the underlying semi-Markov process (x1)t>o- In particular, /(0) is the state 
whose sojourn time distribution has heavy tails. Therefore it follows from Theorem 4.1 of [5] 
that 

C- lim C- lim T l - H \ * Y N ' T \ = [a f g(s u )dB^\ (3.39) 

1 y/W) J < 4 <! Wo J <K1 

for some a > because fi ^ f(0). 

(ii) Let us now define a family of stochastic processes Q^' 7 " = {Q t ' T }o<t<i by 

ft rpl-H rpl/2-H 

Qt = ^{su)Qu du+—==Y t + 7=7= X t . 



* T-.l-.ff rpl/2-H 

1 Y N,T + 

VW) 1 s/W)' 

Since the rate functions are bounded and H > \ 

rpl/2-H 

lim sup sup — X? = 

r-*oo jv o<i<i 7L(T) 

almost surely, and the continuous mapping theorem along with (i) yields 

C- lim £- lim {Qf ,T }o<t<i = {^i}o<t<i 

T— >oo N— >oo 



(iii) Let us put 



Qt = T /TTTtK ^ t -St)- 



Up to a term of the order *^ we obtain 

jVT _ /qE1i{a(^,^ t ) -A(x^,, m )}^ T i-tt ^ T i/2-tt ^ 
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Using the same arguments as in the proof of Theorem 13.131 we thus see that 

lim sup \Q^ ,T — Qt' T \ = i n probability 

N^oo 0<t<l 

for all TsN. Hence the assertion follows from (ii). 

□ 

Remark 3.20 In the case of Markov switching, i.e., when the process xt is a Markov process, we 
obtain standard Ornstein-Uhlenbeck process, i.e., we have that 

C- lim C- lim VT t== ( Sf' T - s t \ = {Zt}o<t<i, 

where Z denotes unique solution to the stochastic differential equation 

dZ t = X'(s t )Z t dt + ag(s t )dB t , 
with B a standard Brownian motion. 

4 Outlook & Conclusion 

We briefly outline two possible avenues of future research: microstructure models of fractional 
volatility and strategic interactions between "big players." 

4.1 Fractional Volatility 

In this article, we suggested a microeconomic approach to financial price fluctuations that is capable 
of explaining the decay of the Hurst coefficient of the S&P 500 index in the late 1990s. We note 
that the evidence of long memory in stock price returns is mixed, there are several papers in the 
empirical finance literature providing evidence for the existence of long memory, yet there are 
several other papers that contradict these empirical findings; see e.g. [8] for an exposition of this 
debate and references. However, long memory is a well accepted feature in volatility (squared 
and absolute returns) and trading volume (see e.g. [23] and |28j). We are now going to illustrate 
how the mathematical results of this paper might also be seen as an intermediate step towards 
a microstructural foundation for this phenomenon. To ease notational complexity and to avoid 
unnecessary technicalities we restrict ourselves to the simplest case where the order rates do not 
depend on asset prices. Specifically, we assume that (after taking the iV-limit) the dynamics of the 
asset price process can be described by a stochastic equation of the from 

where the Gaussian process Y T defined in (|3TTT) converges in distribution to a fractional Brownian 
motion process. In view of the strong approximation of Poisson processes by Browninan motion, 
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and because the rate functions are bounded, the evolution of prices can be described in terms of an 
ordinary differential equation in a random environment generated by a fractional Brownian motion: 



C- lim {Sf }o<t<l = {s*}t<o<l where ds t = X(B^) dt. 
The fluctuations around this first order approximation satisfy 

VT (s? - J* X(Yl) du^j = B + (J* A + (Y„ T ) dv)j - B_ A_(F u r ) du^j 



up to a term of the order ■ Convergence of the Gaussian process Y 1 to fractional Brownian 
motion along with continuity of the rate functions yields 



0<t<l 

Thus, for large T, logarithmic asset prices satisfy 

Sf " jT du + ^[ ^f X +( Y u T ) dB+ --j=f Q ^-(YJ) dB~ 

S jf* A(flf ) du + -^J* ^\ + (BS) dB+ ^-{BS) dB~, 

i.e., the volatility is driven by a fractional Brownian motion process which is independent of the 
Wiener processes B + and B~ . We will further elaborate on the microstructure of fractional volatil- 
ity in a separate paper. 

4.2 Strategic Interactions 

Together with the price taking small investors, it is also possible to incorporate the effects of large 
investors who influence the price. The existence of large agent price effects has been empirically 
described in several papers: [59], [38] and [19] describe the impacts of institutional trades on stock 
prices. In the presence of large agents there is limited liquidity in the market since the holdings 
of the stocks is concentrated in the hands of a few big traders. Trades of "big player's" also affect 
stock prices due to large order sizes. 

4.2.1 Stochastic Equations in Strategically Controlled Environments 

Horst [51] , [52] provides a mathematical framework for analyzing linear stochastic difference equa- 
tion of the form (|2.3p when the dynamics of the random environment is simultaneously controlled 
by the actions of strategically interacting agents playing a discounted stochastic game with com- 
plete information. In [51] we considered a simple microstructure models where small investors 
choose their current benchmarks in reaction to the actions taken by some "big players" . One may, 
for example, think of a central bank that tries to keep the "mood of the market" from becoming 
too optimistic and, if necessary, warns the market participants of emerging bubbles. One may also 
think of financial experts whose recommendations tempt the agents into buying or selling the stock. 
These market participants influence the stock price process through their impact on the behavior 
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of small investors, but without actively trading the stock themselves. It seems natural to assume 
that the big players anticipate the feedback effect their actions have on the evolution of stock 
prices and thus interact in a strategic manner. Under a weak interaction condition, the resulting 
stochastic game has a homogenous Nash equilibrium in Markovian strategies. It turns out that 
the main qualitative feature of the models studied in [40J, [39] and [53], namely namely asymptotic 
stability of stock prices can be preserved even in a model of strategic interactions. However, the 
long run distribution of stock prices depends on the equilibrium strategy and is thus not necessarily 
uniquely determined. Hence, the presence of strategically interacting market participants can be 
an additional source of uncertainty. 

4.2.2 Stochastic Games in a Non- Markovian Setting 

Bayraktar and Poor [7] considered the strategic interaction of large investors and found an equi- 
librium stock price taking into account that the feedback effects of the large investors on the price. 
The large traders find themselves in a random environment due to the trades of small (i.e. price 
taking) investors. In [7j, the institutional investors strategically interact through the controls they 
exert on the coefficients of a stochastic differential equation driven by a fractional Brownian motion. 
Here, the fractional Brownian motion models the effect of the price taking investors on the price. 
It can be argued that the observed stock price is the Nash-equilibrium price that arises as a result 
of the strategic interaction of the institutional investors this random environment. Bayraktar and 
Poor carries out an analysis of stochastic differential games in a non-Markov environment using the 
stochastic analysis for fractional Brownian motion developed in [31]. This analysis can be viewed 
as a first step toward incorporating the feedback effects of the large investors and the strategic 
interaction into the description of the stock price dynamics. 
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